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(page 10) is an approximation which is good only for small 
"exact" equation for a sinusoidal eccentricity is 
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When this equation is used to obtain the "tangent" stiffness of a curved strip, 
the following is obtained in place of equation (All) 
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Use of this equation in constructing figures 2 to 4 leads to results which differ 
from those obtained with the use of equation (All) by a negligible amount. 
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CORRELATION OF THE BUCKLING STRENGTH OF 
PRESSURIZED CYLINDERS IN COMPRESSION OR BENDING 
WITH STRUCTURAL PARAMETERS 
By James P. Peterson 


SUMMARY 


The data on nonpressurized cylinders in Lending of NACA TN 3735 
and the data on pressurized cylinders in compression and bending of 
NASA TN D- 36 O are correlated with structural parameters by using small- 
deflection buckling theory and reduced values for the extensional stiff- 
ness of the cylinder wall. The correlating procedure should prove use- 
ful in future shell-buckling investigations by reducing the number of 
tests required. 


INTRODUCTION 


The current demand for information regarding the buckling strength 
of pressurized shells has prompted several recent test programs. (See 
refs. 1 to 7 .) So far the effort has consisted mainly of obtaining test 
data with little success in interpreting the results of the tests or in 
correlating the data with known structural parameters. The amount of 
test data required in such investigations is large and may be consid- 
erably reduced if a means for correlating data with shell geometry can 
be found. A procedure is developed herein whereby the results of some 
recent cylinder tests are correlated with suitable structural parameters 
by using small deflection buckling theory and reduced values for the 
geometric extensional stiffnesses of the cylinder wall. The reduced 
stiffness values used are devised to provide a basis for correlation 
where geometric deviations from a perfect cylinder at loss of stability 
are responsible for the reduced stiffnesses but may also apply where 
other factors influence the buckling stress. 
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Young's modulus, ksi 

extensional stiffnesses of cylinder wall in axial and cir- 
cumferential directions, respectively, kips/in. 

shear modulus, ksi 

shear stiffness of cylinder wall, kips/in. 
ring spacing, in. 

moment in plate per inch of width of plate, kips 
load in plate per inch of width of plate, kips/in. 
internal pressure, ksi 
radius of cylinder, in. 

length measured along arc of imperfection, in. (fig. 5) 
thickness of cylinder wall, in. 
effective thickness of cylinder wall, in. 


amplitude of imperfection subjected to a load N, in. 
initial value of amplitude of imperfection, in. 


curvature parameter, 
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displacement of load, in. 
arc length of imperfection, in. 
Poisson's ratio 
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|u' x ,M-' y Poisson's ratios associated with extension of the cylinder 
wall in the axial and circumferential directions, 
respectively 

<j normal stress, ksi 

<? cr buckling stress, ksi 


BUCKLING BEHAVIOR OF COMPRESSION CYLINDERS 


A thin-wall cylinder subjected to axial compression usually buckles 
at a fraction of the predicted buckling load for the cylinder. This 
discrepancy between theoretical results and test results appears to be 
largely a function of the radius-thickness ratio of the cylinder and to 
a lesser extent a function of the structural material, the type of con- 
struction, the quality of workmanship, the method of testing, and the 
size of the test cylinder. The reason .for the discrepancy is often 
attributed to local geometric deviations from the perfect cylinder (so- 
called imperfections or eccentricities) and to the associated local 
stress variations which have the effect of rounding off the familiar 
peak in the theoretical load- shortening curve for a perfect cylinder. 

The load that a cylinder with imperfections can carry is considerably 
less than that of a perfect cylinder because the strength of a cylinder 
depends to a considerable extent on the extensional stiffness of the 
cylinder wall. Geometric imperfections can cause an apparent reduction 
in this stiffness because the stretching forces simply tend to increase 
or decrease the imperfections rather than compress or stretch the 
fibers of the middle surface of the cylinder wall. 

The shape of the cylinder just prior to buckling results to a 
great extent from the application of the load in the presence of the 
initial imperfections and, thus, can be expected to resemble the shape 
of the initial imperfections. Therefore, imperfections resembling the 
final buckle configuration would be expected to have the greatest del- 
eterious effect since, otherwise, buckling must be accompanied by a 
change in shape. This phenomenon presumably explains why the discrep- 
ancy between theoretical results and test results for cylinders sub- 
jected to torsion or circumferential compression is relatively small. 

The buckles in the natural buckling mode for such cylinders are very 
long; the buckles are long in the longitudinal direction in the case 
of circumferential compression, and they are long in an oblique direc- 
tion in the case of torsion. The likelihood of imperfections of these 
shapes is not great. On the other hand the theoretical small-deflection 
buckling load of a cylinder in axial compression is very insensitive to 
buckle shape; the buckling mode can consist of many small buckles or a 
fewer number of large buckles. The buckles may be long in the 
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longitudinal direction or long in the circumferential direction. This 
characteristic implies that, in compression, imperfections of a variety 
of shapes and sizes are potentially detrimental. 

The buckling behavior of axially compressed cylinders which buckle 
in the axisymmetric mode appears to be more adequately explained by 
theory than the behavior of cylinders which buckle in the asymmetric 
mode. The axisymmetric mode may be obtained by employing cylinders of 
small radius-thickness ratio or by employing cylinders of larger radius- 
thickness ratio and internal pressurization. For cylinders with small 
radius- thickness ratios, it is sometimes difficult to determine whether 1 

good agreement with elastic theory is obtained or not; cylinders with I 

radius- thickness ratios small enough to insure axisymmetric buckling I 

are so stable that the buckling stress is at least close to the plastic [ 

region of the cylinder material. For cylinders with larger radius- l 

thickness ratios and pressurization, however, axisymmetric buckling can 
be obtained at very low stress levels; in these instances agreement 
between theoretical and experimental results may be satisfactory. (See 
ref. 8.) 


CORRELATING PROCEDURE 


The small- deflect ion buckling load of a cylinder in axial compres- 
sion is strongly influenced by the extensional stiffness of the cylinder 
wall. This characteristic is exemplified in figure 1 which was con- 
structed from computations based on the stability equation (A4) of ref- 
erence 9- The curve labeled G^^Gt in figure 1 was obtained with the 

assumption that all stiffnesses were held constant at their normal value 
except G X y which was varied from the normal value Gt to a very small 

value. The other curve was obtained by three similar but separate cal- 
culations - one in which only E x was varied, one in which only Ey 

was varied, and one in which G^y, E x , and Ey were varied simul- 
taneously by the same amount. The computations indicate that the result 
of a reduced value for any one of the stiffnesses G^y, E x , or Ey is 

a reduced buckling load. Furthermore, a reduced value of one of the 
stiffnesses has about the same deleterious effect on buckling load as a 
reduced value for all three stiffnesses. Inasmuch as the values for 
buckling load shown are not independent of Poisson's ratios associated 
with stretching of the cylinder wall, the values used in the computa- 
tions for figure 1 are given in table I. 
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Development of Procedure 

The phenomena depicted in figure 1 will he used as the basis of 
developing a procedure for correlating the data of references 7 and. 10 
with structural parameters and small deflection theory. For this pur- 
pose the local stiffnesses E x , E y , and G-^ are assumed to be given by 

Ex = %_ Qxy _ le /-j\ 
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where 
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is given by equations (A5) and (All) as developed in the 


appendix for a curved strip. Equation ( A5) applies when the internal 
pressure is zero and may be rewritten for convenience as follows: 


te _ 1 

1 + 6|3 S. 

T> 


( 2 ) 


Equation (All) applies when the internal pressure is not zero and is 
given by 



(3) 


where 



(4) 


and the load N of equation (All) is the hoop tensile force of pres- 
surization pR. The indicated form of the parameters 3 and y was 
chosen to give good agreement between the correlating procedure and 
available test data for constant numerical values of the parameters. 
Additional data may indicate a better form. 


The sweeping assumptions inferred by equations (l) to (3) can be 
justified only by the end result; that is, on the correlation achieved 
and not on the merit of the various assumptions involved. However, 
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some discussion of the assumptions seems warranted. The assumption that 
E x , Ey, and Gy-y can he expressed by a single equation (eq. (l)) is 

made for convenience since figure 1 indicates that cylinder strength is 
relatively independent of whether a single stiffness is reduced or 
whether all three stiffnesses are reduced simultaneously. With this 
assumption, the local buckling stress of compression cylinders is given 
by the familiar stability equation for isotropic cylinders in compres- 
sion with the curvature parameter Z in the stability equation multi- 


plied by 



For the case of moderately long compression cylinders, 


the local buckling stress becomes 



The use of equations (2) and (3) to represent the local extensional 
stiffness of the wall of a compression cylinder just prior to loss of 
stability implies that the major effect of geometric imperfections is 
taken into account by substituting for the small-deflection strain- 
displacement equations, equations which account for the fact that dis- 
placement can occur that does not involve straining of the middle sur- 
face of the shell if imperfections are present. The equations for a 
curved strip are used in place of the corresponding equations for a 
three-dimensional bulge for simplicity; the essential character of the 
resulting equations should be somewhat the same. The actual magnitude 
of the stiffness (that is, the numerical value of 3 and y at loss 
of stability) is established empirically in the following section of 
this paper . 


Application of Procedure 

Nonpressurized cylinders in bending .- The bending test results of 
reference 10 and the nonpressurized bending test results of reference 7 
are compared with equation (5) in figure 2. The two sets of data were 
obtained on similar specimens under similar conditions so the same value 
of 3 should apply for either set. The lower limit of the data is pre- 
dicted fairly well by using 3 = 0 . 00065 ; the average of the data is 
given by 3 = 0.00030. In either case a single value of 3 appears to 
apply over the entire range of data. 


Pressurized cylinders in compression or bending .- The data of ref- 
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erence 7 are presented in figures 3 and 4 as 


plotted against 
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Calculated curves based on equation (5) are also given. The curves were 
drawn by using the numerical values for the parameters 0 and 7 as 
follows : 


Type test 

0 

7 

Compression 

O.OOO65 

1.2 

Bending 

.00030 

3.6 


The value of 0, for the bending tests was taken from figure 2 as being 
representative of the middle of the scatter band for those cylinders. 

A correspondingly reliable value of 0 is not available for the com- 
pression cylinder tests because the number of nonpressurized tests on 
the type of specimens reported in references 7 and 10 is few. The 
value 0 = O.OOO65 for the compression tests was chosen to yield 
buckling stresses which are in reasonable agreement with the data of 
figure 5 and, at the same time, conform with commonly observed ratios 
of compression buckling stress to bending buckling stress. 

The values 7 = 1.2 and 7 = 3.6 were chosen to give good cor- 
relation between the procedure and the middle of the scatter band of 
» the data of figures 3 and 4 at intermediate values of internal pressure. 

The data on cylinders with small values of the curvature parameter Z 
were given less weight in making the choice than data on cylinders with 
larger values of Z. The buckling data on cylinders in bending at high 
values of internal pressure were disregarded in choosing a value of 7; 
however, the data corresponding to strain reversal were considered. 

(See refs. 7 and. 8.) 

Figures 3 and 4 indicate that the stabilizing effect of internal 
pressure is predicted rather well by the correlating procedure; and that, 
at least for the limited data used, reasonable agreement between the 
results of the correlating procedure and experiment can be obtained with 
the use of a single value of 7 for a test series. A correlation of 
the data of references 1 to 6 is not attempted herein. These data evi- 
dently have influencing factors, such as plasticity, type of test, and 
so forth, that did not influence the tests of reference 7* (See dis- 
cussion in ref. 7*) These factors should be evaluated before attempting 
to correlate the results of these tests. 


Discussion 

The correlation achieved in figures 3 and 4 suggests that initial 
imperfections play an important role in the behavior of cylinders at 
low values of internal pressure; whereas at high values of internal 
pressure, any detrimental imperfections that might have existed prior 
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to pressurization are pulled out by the pressurization stresses so that 
they play a lesser role. The pressurization stresses themselves create 
radial prebuckling deformations or imperfections. These deformations 
are apparently such that they are not very detrimental as far as initi- 
ating asymmetric buckling is concerned. Otherwise pressurization might 
have a destabilizing effect rather than a stabilizing effect. The pre- 
buckling deformations caused by pressurization are such that the effec- 
tive extensional stiffness of the cylinder wall in the longitudinal 
direction should be decreased. For this case asymmetric buckles ordi- 
narily would tend to be long in the axial direction and short in the 
circumferential direction. The prebuckling deformations are the oppo- 
site of this; that is, they are short in the axial direction and long 
in the circumferential direction, and perhaps not very influential in 
initiating asymmetric buckling. This argument would seem to be strength- 
ened by the fact that the asymmetric buckles which finally appear when 
a highly pressurized cylinder is loaded to large values of shortening 
have about the same axial wave length as the axisymmetric prebuckling 
deformations which they replace. 

The correlation achieved in figures 2 to 4 indicates that changes 
in buckling coefficient with changes in radius- thickness ratio of non- 
pressurized cylinders as well as the stabilizing effect of pressuriza- 
tion in pressurized cylinders is predicted rather well by the correlating 
procedure. The procedure should therefore prove useful in eliminating 
the need for much of the experimental data that are ordinarily required 
in the establishment of design curves for a new structural material, 
for a new method of construction, or for a different method of loading. 

It can be particularly useful in cases where the effect of pressure 
stabilization is desired on cylinders having proportions such that a 
single cylinder cannot be used for obtaining data at more than one 
value of internal pressure. 

CONCLUDING REMARKS 

A procedure is presented whereby the results of recent cylinder 
tests are correlated with structural parameters by using small deflec- 
tion buckling theory and reduced values for the extensional stiffness 
of the cylinder wall. Values of the reduced stiffness were devised on 
the basis that geometric imperfections may cause substantial reductions 
in the effective extensional stiffness of the cylinder wall which, in 
turn, may cause substantial reductions in cylinder strength. The cor- 
relating procedure may be useful in future experimental investigations 
by eliminating the need for much of the data that is ordinarily required 
to establish empirical design curves. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., July 19, 19^0. 
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APPENDIX 


DETERMINATION OF THE STIFFNESS OF A CURVED STRIP 


The stiffness of a strip as shown in figure 5 with an eccentric- 
ity w Q and subjected to the action of the loads N is desired. It 
will be obtained by Castigliano' s theorem which states that the dis- 
placement of the load A is related to the strain energy of the strip 
according to the formula 



The internal bending moments M and middle-plane stresses a required 
for equilibrium are 

M = Ny (A2) 

and 

a = — cos 0 (A3) 

The symbols y and 0 are defined by figure 5* 

When equations (A2) and (A3) are substituted into equation (Al) 
and the indicated differentiation performed, the stiffness Et e can 

be written after some simplification as 



The two terms in the denominator of the right side of equation (A4) can 
be approximated reasonably well with rather simple expressions. The 
first term is essentially unity for shallow eccentricities of the type 
of interest here, and the second term will be recognized as the ratio 
of the moment of inertia of the curved strip about the line x-x 
(fig. 5 ) to the moment of inertia of a flat strip. This term can be 
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approximated closely by 6^p.J • With the use of these simplifications, 

the effective stiffness can be written in terms of the eccentricity w Q 
as 



(A5) 


Equation (A5) gives the initial stiffness of the strip. That is, it 
gives the stiffness of the strip as the load N approaches zero, hence 
the subscript o on w. For this case (N ^»0), the stiffness depends 
mainly (according to eq. (A5) depends completely) on the depth of eccen- 
tricity w 0 and is relatively independent of the length or shape of 

the eccentricity. Actually the stiffness is not completely independent 
of these characteristics but equation (A5) neglects these effects. 


In order to determine the stiffness of the strip (fig. 5) as a 
function of load, further considerations are necessary. As the load 
N is increased, the eccentricity w Q will be reduced and the stiff- 
ness Et e will be increased. For very large loads the stiffness Et { 
will approach Et. The equation expressing the behavior of the strip 
subjected to an incremental load dW is 0 


dA = dN 


Et* 


= dN ii- JL = 
Et to 


dN 


Et 


1 + 6 


(?) ! 


-j^*(A6) 


where w is the instantaneous value of eccentricity depth, not the 
initial value. The displacement A is the result of a geometric change 
in shape of the eccentricity and a stretching of the fibers of the strip 
which can be expressed in equation form, by assuming a sinusoidal eccen- 
tricity, as 




A - i ^ 2 ■ “ 2 ) + e (A7) 

The first term of the right side of equation (A7) gives the displace- 
ment resulting from a geometric change in shape of the eccentricity 
from one with an amplitude of w 0 to a more shallow one with an ampli- 
tude of w. The last term gives the displacement resulting from a 
stretching of the middle surface fibers of the eccentricity. Equa- 
tions (A7) can be rewritten as 
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V * J * 2 t 2 «2 Et 3 


(A8) 


The substitution of equation (A8) into equation (A6) yields the desired 
differential equation 


I 





1 + 



24A A + 2hNA 2 
rt 2 t 2 * 2 Et3 


(A9) 


The solution of equation (A9) for the boundary condition A = 0 when 
N = 0 is 


n 2 t 2 
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n 2 Et5 


(A10) 


A comparison of equation (A10) with equation (A8) 

2hNA 2 

^2 ' 
e 


yields 
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or 




(All) 


Equations (A5) and (All) are used in the body of this paper to 
affect a correlation between data and theory for nonpressurized and 
pressurized cylinders, respectively. 
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Figure 2.- Application of correlating procedure to nonpressurized 

bending tests. 
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Figure 3.- Concluded 
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(a) Cylinder 5. 



(b) Cylinder 6A. 

Figure 4.- Application of correlating procedure to pressurized bending 
cylinders of reference 7* Symbols with tails denote strain-reversal 
stresses . 
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(c) Cylinder 7A. 
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(d) Cylinder 8. 
Figure 4.- Continued. 
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(b) Differential element. 

Figure 5*- Dimensions of curved strip. 
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